Quasi-particle approach for general lattice Hamiltonians 
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In many condensed-matter systems, it is very useful to introduce a quasi-particle approach, which 
is based on some sort of linearization around a suitable background state. In order to be a systematic 
and controlled approximation, this linearization should be justified by an expansion into powers of 
some small control parameter. Here, we present a method for general lattice Hamiltonians with 
large coordination numbers Z ^ 1, which is based on an expansion into powers of 1/Z. In order to 
demonstrate the generality of our method, we apply it to various spin systems, as well as the Bose 
and Fermi Hubbard model. 
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I. INTRODUCTION 

Apart from a few exceptions (mostly in one spatial 
dimension), the majority of Hamiltonians describing in- 
teractions in condensed matter and other fields is not ex- 
actly solvable. Thus, it is often very useful to introduce 
a suitable quasi-particle approach - i.e., to consider lin- 
earized excitations (quasi-particles) around a given back- 
ground (e.g., the ground state). However, in order to 
avoid neglecting contributions which are of equal mag- 
nitude as the terms one keeps (or even larger), such a 
linearization should be based on the expansion into pow- 
ers of some small control parameter. 

As an example, Bose-Einstein condensates of dilute 
atomic gases are often described in terms of a mean-field 
expansion where the full field operator is split up via 



(1) 



into the condensate wave- function tpcitif) and the re- 
maining (linearized) quantum or thermal fiuctuations 
S^{t,r). The above split can be justified if a large 
number ^ 1 of atoms share the same wave-function 
^pc{t,r), since it scales with \/N in this case [H, HI- 

As another example, spin systems are often described 
using the large-S" expansion (spin- wave theory) [l3| . 
This approximation can be motivated via the Holstein- 
Primakoff transformation which maps the spin operators 
to bosonic creation and annihilation operators a} and a 



S+ = aW2S-a'ia = a^V2S+ 0{l/\^) (2) 

and S- = S\, as well a.s Sz = S — at a. For large spins 
S, the system is then approximated by non-interacting 
bosons. The reliability of the l/S'-expansion for 5* = 1/2 
systems is another question, see Section [IVI below. 

In the following, we present a general formalism for 
the introduction of quasi-particles for lattice Hamiltoni- 
ans which is based on an expansion in powers of 1/Z, 
where Z is the coordination number of the lattice. The 



paper is organized as follows. Section HIl provides the key 
elements of the 1/Z method. In section lHll we derive the 
formalism for quasi-particle excitation. We illustrate the 
usefulness of these concepts for calculating the excitation 
spectrum and the phase transition boundaries for various 
known model descriptions: the Hcisenberg model in sec- 
tion |lVl the Ising model in section |Vl the Bosc-Hubbard 
model in section I VII and the Fermi-Hubbard in section 
IVIII The last section IVHIl is devoted to the conclusions. 



II. BASIC FORMALISM 

We consider a general lattice Hamiltonian 

^=|E^a.- + E^a.' (3) 

consisting of on-site terms iJ^ and two-site coupling or 
tunneling terms Hfiu (with H^^ = 0) where ji and v label 
the lattices sites. For simplicity, we assume that the 
have all the same form (and similarly the H^u), i.e., we 
have (discrete) translational and rotational invariance of 
the underlying lattice. Finally, the coordination number 
Z counts the number of tunneling or interaction neigh- 
bors /X of a given lattice site v and we consider the limit 
Z ^ 1. The time-evolution of the density matrix p of 
the total lattice is given by the von Neumann equation 

1 
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(4) 



with the Liouville super-operators C^^p = [H^^,p\ and 

C^p = [H^,p]. Now we introduce the reduced density 
matrices p^ (or p^ti,) of a given lattice site p, (or p, v) by 
tracing over all sites except p (or p, v) 



/5p = Trj/{/5} , p^^ = Tvyi^/{p} , etc. 



(5) 



Note that p^j is a linear operator acting on the Hilbcrt 
space T-L^i of one lattice site p while p^^^, is a linear oper- 
ator acting on 7^^ ® T-Lu etc. Now we may split up those 
matrices into correlated and uncorrelated parts 



' lralf.schuetzhold@uni-due.del 



Pl^v 



Pp.p^ ■ 



(6) 



2 



and p^^A - P^^A + P^iv P\ + P,xA Pv + Pi/A + Pt^Pi^P^ 
etc. From Eq. (U), we get the equation of motion for the 
reduced density matrix 

^dtP, = I E Ti-. { C.%{PIT + M] + , (7) 

where C^^^ = + 

The above equation is stih exact - but in order to solve 
it and to obtain p^, we would need the correlations /5^°"'. 
Similarly, the equation for the two-point correlations ff°" 
contains the three-point correlator /5™" and so on. How- 
ever, for large coordination numbers Z ^ 1, one can 
solve this set of equations approximately and iteratively. 
If we start in a separable (i.e., uncorrelatcd) initial state, 
the correlations satisfy - at least for a finite period of 
time - the following hierarchy 0, Q : 

p,= (ZO) , = O {1/Z) , pI°J: = O (1/Z2) (8) 

and so on. This allows us to obtain the approximate 
equation for the on-site density matrix 

i9tP^ = I E { C^^Jt^Pn} + C,,p,, + 0{1/Z) (9) 

The zeroth-order solution p^ of this equation corresponds 
to the local (e.g., Gutzwiller) mean-field ansatz. Insert- 
ing this solution and performing an analogous procedure 
for two sites, we get 

idtplT = £,PlT + \ £,upIpI - f Tr^ { CIjIpI] 

+i { ^UpITpI + pITpD} 

K/ 

+{p^v)+ 0{l/Z^). (10) 

This equation determines the quantum fluctuations 
around the background solution pj^ and will be the start- 
ing point for our quasi-particle picture. 

III. FACTORIZATION 

As it turns out, Eq. (jlOp can be simplified considerably 
by effectively factorizing it. To this end, let us introduce 
the propagator W^^(t,io) from an initial lattice site v at 
time to to a final lattice site p at time t which corresponds 
to the mapping 

Mh) ^ i.(t) = E TrM{w^;(i>^o)iM(io)} .(11) 

Here W'^(i,to) is a linear operator acting on ® Hi, 
where and Hi, are two different Hilbert spaces (initial 
and final lattice sites) even for p, = v. 

To achieve consistency for t = to-, wc impose the initial 
condition 

W;{t, t) m^) = 5^-^ , (12) 



where the |n^) form a complete basis of H^ and simi- 
larly the \na) for Ha- If we now postulate the following 
effective equation of motion for the propagator 

+1 E Ti-K { [plW^it, to) + W^{t, to)pi) } (13) 

K. 

and insert it into the factorization ansatz 

PlTit) = E Tr.^ {w^{t,to)Wl'{t,to)plJ{to)} + 

al3 

t 

where Qq/j denotes the sources term in Eq. ([TO)) 

= * f Tr. { £f,p°p° } + ^ | Tr, { ff,p°pO } 

~^C%plfPp, (15) 

we find that we obtain the solution of Eq. (fTO)) to first or- 
der in 1/Z. Ergo, instead of solving Eq. ((TU)) directly, we 
may solve the simpler equation (|13p instead and recover 
the correlations p'^^"{t) as a bilinear combination 

Interestingly, the effective equation ([T3l) for fixed a 
(the initial lattice site) is formally equivalent to a direct 
linearization of the on-site equation To see that, let 
us set Pfj, = p^ + Spfj, and linearize around the stationary 
solution pj^. The resulting equation 

idt6p^ = C^5p^ + I E (plSp^ + Sp^pl)} 

K 

(16) 

is formally equivalent to p3l) . Accordingly, the solution 
of this equation is given by (jlip 

Sp,,{t) = Tr. [w;:{t,to)dMh)} , 

which confirms that W'^i^.to) is the propagator. 

IV. HEISENBERG MODEL 

The Heisenberg model is the first basic lattice model 
to which the 1/Z approach applies. The Hamiltonian is 

H^-^Y.^,.S^-S,. (17) 

Here S^i = {^^: ^fn a-^e spin operators at the lattice 
site p satisfying the usual SU (2) commutation relations 
[S'^,S'^] = i5^„e°'^''S^. J denotes the coupling strength. 
The lattice structure is encoded in the adjacency matrix 
T^i, which equals unity if p and f are tunnelling neigh- 
bours (i.e., if a particle hops from p to v) and zero other- 
wise. For simplicity, we restrict ourselves to spin 5* = 1/2 
systems, but our results can be generalized easily. 
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A. Ferromagnetic case J > 

The Heisenberg model (|17p is invariant under internal 
(spin) rotations, but its ground state breaks this sym- 
metry spontaneously, i.e., we have a set of degenerate 
ground states. For positive J, any state where all the 
spins point into the same direction is an exact ground 
state. For simplicity, let us take the spin-up state: 



pI = lt>,(tl 



(18) 



Introducing the time-propagated spin operator 



S.{t) = ^ Tr^ 



[w;{tM)S^{h)] , (19) 



according to ([TT|) . we obtain the following closed opera- 
torial equation of motion: 



Z 



(20) 
(21) 

where we have used the notation — Sf, i: iSf,- As- 

M M _ 

suming (discrete) translational symmetry of our lattice, 
we may Fourier transform these equations 



(w T J[l - Tk] 



= 0, 



(22) 



where Tk is the Fourier transform of the adjacency matrix 
Tf^u with the normalization Tk=o = 1- From the above 
equation, we may directly read off the spectrum 

Wk = ± J(l - Tk) . (23) 

As expected, the magnons in the ferromagnetic phase are 
gapless (Goldstone modes) and their spectrum scales as 
Wk ~ for small k. Note that the same spectrum can 
be derived via the l/S-expansion [l3l | mentioned in the 
Introduction. However, for 5* = 1/2, the 1/Z-cxpansion 
appears better justified than the 1 / S'-expansion. 

The positive frequency associated to the excitation op- 
erator S^^ means that it costs energy to flip one spin up. 

On the contrary, has a negative energy which means 
that adding a spin up to a spin polarized magnet liber- 
ates energy. This process is not possible because already 
all the spin are up. Therefore, such a release of energy 
is possible provided excitations are created in pair: a 
creation of positive energy followed by a creation of neg- 
ative energy so that the net energy is positive. A switch 
in the sign of J changes the sign of the frequencies and 
announces the change to the antiferromagnetism phase. 
In this case, the system becomes instable because energy 
is released by changing a spin up into a spin down. 



B. Anti-ferromagnetic case J < 

For negative J = — Ja, the ground state(s) cannot 
be determined exactly in general. Nevertheless, under 



certain conditions, we can get a good approximation. In 
order to avoid problems with frustration, we assume that 
we have a bi-partite lattice, i.e., that we can divide the 
lattice into two sublattices A and B such that for each 
site /I G .4 all the adjacent neighbouring sites v belong to 
B and vice versa. Then, for large Z, we may approximate 
the anti-ferromagnetic ground state by the Nccl state [l^ 



(24) 



A similar reasoning as in the ferromagnetic case leads to 
the following operatorial equations 



(25) 



{^^,±JA)S^^TY.'^^St if ^i^A, (26) 
{idtTJA)S^ = ±Y.'-^^St if f^eB. (27) 



After a Fourier transform, we deduce the following spec- 
trum for the magnons in the anti-ferromagnetic case 



(28) 



As expected, the spectrum is still gapless but scales as 
Wk ^ |k| for small k. Again, the same spectrum can be 
derived via the l/S-expansion [l^ 

Also the frequencies go to zero at the 
antiferromagnetism-ferromagnetism transition. At 
the difference of ferromagnetism, it costs energy to 
change the spin up from spin down in sublattice A 
and the other way around for sublattice B. But this 
process reverses once the frequencies go to zero at the 
antiferromagnetism- ferromagnetism transition. 



V. QUANTUM ISING MODEL 

As our second example, let us inmvestigate the quan- 
tum Ising model in a tranverse field with the Hamiltonian 



z 



(29) 



This model displays a quantum transition from the para- 
magnetic phase (where the transversal magnetic field B 
dominates) to the ferromagnetic or anti-ferromagnetic 
phase, which breaks the Z2 spin-flip symmetry. 

To lowest order in 1/Z, we employ the following vari- 
ational ansatz for the ground state 

^0 = E [h^p'',] + I E T^m- {^p./5°P°}(30) 

where the reduced density matrix is given by a pure state 
P° = IV'>m(V'I with 1^)^ = ctl t)^ + c^l Da.- The varia- 
tional ground state energy per site reads 



El 

N 



J, 



B 



cth - |c;| )" - -7^{c\cx + c^cD . (31) 
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Minimizing this energy, we find that the ground state al- 
ways has (iS'po = 0, i.e., we can choose the amphtudes 
c-|- and ci to be reaL To lowest order in 1/Z, the crit- 
ical point is at J = B. For J < B, the magnetic field 
controls the orientation of the spin so that the state is 
paramagnetic with (S'j^) = and {S^) = 1/2 (assum- 
ing i? > 0). For J > B, on the other hand, the fcrro- 
magnetism case is associated to a permanent magnetiza- 
tion with a non- vanishing ferromagnetic order parameter 

(%> = ±^i~Byjy2. 

In complete analogy the Heisenberg model, we may 
derive the effective quasi-particle equations 



(32) 



B 



^dtS^^ = -iS+ - S-) . 



This gives the excitation spectrum 



Wk = ± 



^AJ^S~;)l + B^-2B.m{SS) 
For the paramagnetic phase, this simplifies to 



LUi, = ±y^B^ - BJTk. 



(33) 



(34) 



At the difference of the Heisenberg models, here tran- 
sition occurs for J = i? when the frequencies becomes 
imaginary. For comparison, the spectrum in ID obtained 
from exact diagonalisation is given by 



= ±v/B2 - BJTk + JV4 



(35) 



and is well approximated by the lowest order result 
Ea. (|M|) for large B. In ID however transition occurs 
for a zero frequency at B = J/2. Let us note here that, 
although the 1/Z expansion is supposed to work very 
well in the limit of large dimensions, it gives nevertheless 
qualitatively good results for ID. 

In the ferromagnetic case, the spectrum is 



(36) 



and corresponds in the limit of weak B to two kink exci- 
tation 2wi£( of the exact ID spectrum Eg. ((55)) . 



VI. BOSE-HUBBARD MODEL 

After discussing two examples for spin systems, let us 
apply the hierarchy discussed above to the Bose-Hubbard 
model, see also Q • This model describes (identical) 
bosons hopping on a lattice with the tunnelling rate J. If 
two (or more) bosons occupy the same lattice site, they 
repel each other with the interaction energy U. Alto- 
gether, the Hamiltonian reads 

^ = -|E^m-^1'^- + |E"a<(^m-1)- (37) 



Here and bi, are creation and annihilation opera- 
tors at the lattice sites /z and v, respectively, satisfy- 
ing the bosonic commutation relations [6^, bl] — 6^i, and 

WfiiH'] = \bn^^v\ — 0. As before, the lattice structure is 
encoded in the adjacency matrix T^^ which equals unity 
if ^ and v are tunnelling neighbours (i.e., if a particle can 
hop from /z to v) and zero otherwise. Finally, ~ 
is the number operator and we assume the Mott phase 
P% ^ following. 

In complete analogy to the spin operators in the previ- 
ous Sections, we may define effective particle and hole 
operators via /i^ = = |1)^(2|. Time- 

propagating these operators according to ([TT|) . we deduce 
the following operatorial equations (9|: 



Z 



(38) 



Thus, the eigenvalues obtained from this system are 



U - JT^ ± ^JV^ - 6JUTk + J^T^ 



(39) 



They correspond to the ones obtained in the random 
phase approximation 0, [HI . Inspection of the frequen- 
cies shows that they become imaginary for J/U — 3— 2a/2 
indicating the onset of instability of the Mott state to- 
wards a superfluid state. 



VII. FERMI-HUBBARD MODEL 

As our final example, let us consider the Fermi- 
Hubbard model, which is somewhat similar to the Bose- 
Hubbard model, but significantly more complicated. The 
Hamiltonian has a quite similar form as above 



(40) 



but now the cj^ ^ and c,^,s are fermionic creation and an- 
nihilation operators and we have two different fermion 
species labeled by their spin s =t and s =|. Depending 
on the values of J and U as well as the filling factors 
(np and {njj}, there are various phases which all require 
a different ansatz for (fl, . 

Let us consider the case of large and positive U with 
small fillings (np ^ 1 and {njj} ^ 1 and where the 
temperature T satisfies U ^ T ^ J . In this case, double 
occupancy |ti)^ = cj^ |cj^ 4, |0)^ can be neglected and any 
possible spin ordering is washed out thermally. We may 
then use the following ansatz 

P° =Po|0),.(0| +Vi\ i)A I +Ptl t)p(t I , (41) 

where Po + Pi + Pt ~ ^ with (np ~ and (rip = p^. 
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This ansatz provides a solution of the lowest (zeroth) 
order in XjZ and we find that the quasi-particle modes 
decouple into two branches for the two spin degrees of 
freedom. Let us analyze the spin-up excitation (the re- 
sults for the other spin-down branch are analogous) . The 
relevant perturbations part of the density matrix read 

% ^^^A'.tl t>,.(0| + «p.tl TDmU I 

where = c^,t(l - and = c^.t^^",; 
five particle and hole operators respectively. Using this 
ansatz in Eq. for the Hamiltonian (PH]). we obtain 
the following set of equations 

K 

{idt - C/)wM,t = ""^ H T).Kn("«,t + ^^^.t) ■ (43) 

K 

After a Fourier transform, we obtain the following Hub- 
bard expression for the eigen-frequencies 

4 = ^ ± 2 V + " 4C/jrkn • (44) 



VIII. CONCLUSIONS 

Via a controlled expansion into inverse powers of the 
coordination number 1/Z leading to a hierarchy of cor- 
relations [1, Q , we establish a simple and elegant way to 
obtain quasi-particle operators and spectra. We demon- 
strated the general applicability of this method by means 
of several examples such as spin systems, the Bose and 
Fermi-Hubbard model. Even though some of the results 
can also be obtained by other means - such as the 1/S- 
expansion - the 1/Z-cxpansion developed here is in sev- 
eral cases (such as for S = 1/2) better justified than the 
1/5-cxpansion, for example. 
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